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ABSTRACT 

If the stellar population of the bulge contains black holes formed in the final 
core collapse of ordinary stars with M ^ 30M Q , then about 25,000 stellar mass 
black holes should have migrated by dynamical friction into the central parsec 
of the Milky Way, forming a black hole cluster around the central supermassive 
black hole. These black holes can be captured by the central black hole when 
they randomly reach a highly eccentric orbit due to relaxation, either by direct 
capture (when their Newtonian peribothron is less than 4 Schwarzschild radii), 
or after losing orbital energy through gravitational waves. The overall depletion 
timescale is ~ 30 Gyr, so most of the 25,000 black holes remain in the central 
cluster today. The presence of this black hole cluster would have several 
observable consequences. First, the low-mass, old stellar population should 
have been expelled from the region occupied by the black hole cluster due to 
relaxation, implying a core in the profile of solar-mass red giants with a radius 
of ~ 2pc (i.e., 1'). The observed central density cusp (which has a core radius 
of only a few arc seconds) should be composed primarily of young 1 Gyr) 
stars. Second, flares from stars being captured by supermassive black holes in 
other galaxies should be rarer than usually expected because the older stars 
will have been expelled from the central regions by the black hole clusters of 
those galaxies. Third, the young (^ 2 Gyr) stars found at distances ~ 3 — 10 pc 
from the Galactic center should be preferentially on highly eccentric orbits. 
Fourth, if future high-resolution A"-band images reveal sources microlensed by 
the Milky Way's central black hole, then the cluster black holes could give rise 
to secondary ("planet-like") perturbations on the main event. 



Subject headings: black hole physics - Galaxy: center - Galaxy: kinematics and 
dynamics 
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1. 



Introduction 



The measurement of proper motions and radial velocities of stars within the central 
parsec of the Galaxy has led to the conclusion that a black hole of mass 3.0 ± 0.3 x 1O 6 M 
is present in the center (Eckart & Genzel 1997; Genzel et al. 2000). There is also increasing 
evidence that massive black holes are found in the centers of other galaxies (Richstone et 



The central region of the Galaxy is also peculiar because the relaxation time among 
stars can be shorter than the age of the Galaxy owing to the high density. The process of 
relaxation leads to a stellar cusp, which has a density profile p oc r -7 / 4 when all the stars 
have the same mass (Bahcall & Wolf 1976). Several interesting physical processes take place 
among the stars in this cusp: stars can come close enough to physically collide with each 
other, and they can also come sufficiently close to the black hole to be tidally disrupted 
or swallowed (e.g., Frank & Rees 1976; Lightman & Shapiro 1977; Quinlan, Hernquist, & 
Sigurdsson 1995; Sigurdsson & Rees 1997). 

One of the consequences of the relaxation is that the most massive objects will sink to 
the center of the stellar cusp. Among an old stellar population, the most massive objects 
should be black holes formed in the final core collapse of massive stars. We assume in this 
paper that most massive stars with M ^ 30 M Q produce black holes, most with a mass of 
7M (Bailyn et al. 1998). The high mass of these black holes implies that their dynamical 
friction time to move to the center of the Galaxy is shorter than a Hubble time over a much 
larger volume than the one where ordinary stars have a short relaxation time. We will find 
in §2 that this should lead to the formation of a cluster of stellar black holes around the 
central supermassive black hole (hereafter "Sgr A*"), and that other stars are ejected from 
the region occupied by this cluster. In §3 we discuss the rate at which the black holes in 
this cluster are captured by Sgr A*, and we find that most of the black holes should still 
be present in the cluster. Several observable consequences of the presence of this black hole 
cluster are discussed in §4. 



The deprojected light profile in the inner kpc of the Galaxy scales as r -1 ' 8 , while the 
predicted profile around a massive black hole (r ^ 1 pc) scales as r~ 7 / 4 (Bahcall & Wolf 
1976). For simplicity, we therefore adopt a density profile p(r) 



al. 1998). 



2. 



Cluster Formation 



p(r) oc r 



-7/4 



(1) 
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From the model fit of Genzel et al. (2000) to the velocity dispersion data, we find that the 
total mass within r = 1.8 pc is 2M cbh (see their Fig. 17), where 

M cbh = 3 x 10 6 M (2) 

is the mass of Sgr A*. Hence, the total distributed mass inside 1.8 pc is M cbh , and the 
density profile is 

P*( r ) = ITT- — (—) ' r o = L8 P C - ( 3 ) 



167r rg Vr 

We assume that this density profile is entirely composed of stars, brown dwarfs, and 
stellar remnants. To calculate the mass fraction of black holes, r/ hh , we use the following 
initial mass function. For the range M > 1 M , we adopt a Salpeter law dN/dm oc mr a 
with a = 2.35. For 0.7 M < m < 1M , we adopt a = 2 from Zoccali et al. (2000). 
For 0.05 M Q < m < 0.7 M , we adopt a = 1.65 by correcting the result of Zoccali et al. 
(2000) for binaries according to the adjustment of Gould, Bahcall, & Flynn (1997) and by 
extending the power law beyond the last observed point at 0.15 M . We cut off the mass 
function at m ~ 0.05 M in accordance with the prelimary indications from microlensing 
(Han & Gould 1976). We assume that all progenitors with masses 1-8 M have become 
0.6 M Q white dwarfs, those with masses 8-30 M have become 1.4 M neutron stars, and 
those with masses 30-100 M Q have become 7M black holes. We then find 

77 bh = 1.6%, (m) = O.23M , (4) 

where (m) is the mean mass of the population. 

Once this population is formed, the black holes will sink toward the center on the 
dynamical friction timescale (Binney & Tremaine 1987) 

nA AnGpGrn^ r duW/(u) (5) 
v 6 Jo 

where m hii = 7 M Q is the mass of the black hole, v is its velocity, f(u) is the velocity 
distribution of the ambient stars, and the integral gives the fraction of ambient stars with 
speeds below v. For a Gaussian velocity distribution with dispersion a, and a typical black 
hole speed v 2 = 3a 2 , the value of the integral is 0.54. For the Keplerian part of the potential 
(r < r ), InA = ln(M cbh /m bh ) = 13. For r > r , InA rises slightly but we ignore this in the 
interest of simplicity. We evaluate idf,o> the fiducial dynamical friction time at r , 

t df , = 1.4Gyr (6) 

and note its scaling in the two regimes 

/ r \ 1/4 / r \ 17/8 

tdf = *df,o ^— J (r < r ), tdf = ^f.oj^— J (r>r ). (7) 
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Hence, after a time t, all the black holes that were originally within a radius r will collect 
in a cluster near the center, where r is given by, 



r 
r 



4t d f,o 



(t < 4t 



df.OJ 



r 
r 



17/ t 



17/ 
7 U 



df,0 



4 +1 



8/17 



(t > 4t df , ). (8) 



This implies an infall rate of black holes 



dN, 



5 ?7bh M chh ( t 



dt 



7^—) (t < 4tdf,o) 

4tdf ,0 / 



dN hh = lQ Vbh M chh 
dt 7t dli0 m hh 



7 v^df.O 



-7/17 



(t > At, 



df,0j 



(9) 



If we assume that the bulge formed at a time tbuige ~ 10 Gyr, then from equation (J|), 
all the black holes within a radius r d f = 5 pc will have migrated to the center by now. The 
cumulative total and current rate of precipitation are therefore 



A 7 i 



bh 



2.4 x 10 4 , 



bh 



dt 



2.9 Myr" 



(10) 



In other words, provided that our assumption of the fraction of massive stars that form 
black holes in their final core collapse is correct, we must conclude that a large number of 
stellar black holes, with a total mass of ~ 5% of the Sgr A * mass, have migrated to the 
center and, unless they have subsequently been captured by Sgr A *, they should have formed 
a cluster of black holes in the center of the stellar cusp. 

As the black holes precipitate, they start dominating the total density in some central 
region at some point, and then the low-mass stars are expelled from this region over a 
relaxation time. Assuming that most of the energy is lost from the cluster by direct 
capture of black holes near the center, the black holes should also relax to a density profile 
proportional to r -7 / 4 , for which the outward energy flow is constant. As this energy flow 
is transmitted to the low-mass stars outside the black hole cluster, the cluster will need to 
expand and push out the low-mass stars. To derive the relative density of the black hole 
profile compared to the empirically normalized stellar profile, we invoke the steady state 
energy-flow condition between two species of stars A and B of mass and mg, which 
dominate the total density pa and ps at radii and re, respectively. The total energy 
at radius r is proportional to pa 2 r 3 , and the relaxation time is proportional to cr 3 /(pm). 
Therefore, the constant energy flow condition yields 



PA<y A r A 



PBO- 2 B r% 



a 3 /(p A m A ) G- 3 /(p B m B ) 



(11) 
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Making use of the Kepler-potential relation a 2 oc r, this implies, 

pa f^r 7/ V^r 1/2 



Pb \r B J \m B J 



(12) 



Thus, the mass density of black holes in the central region is below that implied by 
extrapolating equation (|3|), by a factor ((m)/m bh ) 1//2 , that is, Pbh( r ) = 0.18p*(r). Hence, 
if all of the black holes precipitated over the lifetime of the galaxy from a radius of 5pc 
remain in the cluster at present, the cluster should extend over a radius rt,h, 

rbh = 7 x 5pc = 0.7pc. (13) 
V (to) / 

The timescale to achieve this expansion is the relaxation time in the expanded (lower 
density) cluster, which is ~ 3 times longer than the dynamical friction timescale because of 
the lower density. 



3. Rate of Capture of the Black Holes 

In the previous section we found that about 24,000 stellar black holes should have 
migrated to the center of the stellar cusp around Sgr A*. We now address the question of 
the rate at which these black holes will be removed from the cluster by coalescing with Sgr 
A*. If this rate is low, most of the black holes should be in the cluster at present. If the 
rate is high enough, then many fewer black holes will be present, and a balance between the 
rate at which black holes are precipitating in the cluster by dynamical friction and the rate 
at which they are being captured should be established. 

The dominant process by which black holes will be eliminated is by a random walk 
into a highly eccentric orbit as their orbits change over the relaxation timescale, from which 
they can be captured by Sgr A*. This process was first studied by Frank & Rees (1976). 
In the case of stars, tidal disruption can eliminate them from the cluster once they come 
close enough to Sgr A*; obviously, orbiting black holes will be eliminated only when they 
are swallowed by Sgr A*, possibly after having lost orbital energy by emitting gravitational 
waves. 

Before describing in more detail the mechanism by which black holes are captured, we 
need to discuss the process of orbital diffusion by which black holes will migrate into the 
eccentric orbits from which they can be captured. 
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3.1. Orbital Diffusion 

A black hole can be captured by Sgr A* from an orbit of any semimajor axis, provided 
that its peribothron q = a(l — e) is sufficiently small. This will lead to a distribution of 
black holes in phase space that is strongly depleted at eccentricities very close to unity, and 
diffusion of black holes will take place toward orbits of decreasing peribothron. In order to 
investigate quantitatively this black hole migration, we first evaluate the diffusion tensor in 
velocity space. We sketch the derivation here and leave the details and the justifications for 
the various approximations to Appendix A. 

The diffusion equation is given by 

where / is the phase-space density, j is the "velocity current density" , and km is the diffusion 
tensor. By symmetry, K ki = diag( ), where k± and n r are the components of k 

perpendicular and parallel to the radial direction. In general, the diffusion tensor depends 
on the spatial position and the velocity. A useful physical interpretation of the diffusion 
tensor components is that, over a small interval of time St, the rms change in the velocity of 
a black hole in a direction perpendicular to its initial velocity is equal to (2k± St) 1 / 2 . The 
total rms change in any direction is therefore (2tr(K)5t) 1 ^ 2 (where tr means the trace), and 
the relaxation time is of order v 2 sc /[2tr(K,)]. 

We assume that the unperturbed (by black-hole capture) phase-space density, / , is 
a function only of the energy, and hence find, for a p oc r~ a density profile in a Kepler 
potential, 

f (u,r)=g(^-)h(r), g{x) = (1 - x) a -^Q{\ - x), (15) 



v 

CSC 



w , = (3-a)a! A bh / r X 3 / 2 '" 

{r> 8n 2 (a - 3/2)1(1/2)! {2GM chh r hh fl 2 \rJ ' 1 ' 

where v esc (r) is the local escape velocity, Abh is the total number of black holes within 
a radius r^, M c bh is the central mass, and G is a step function. Thus, for a = 7/4, 
fo oc r _1 / 4 (l — v 2 /v csc) 1 ^ 4 - In Appendix A, we show that for a = 3/2, the velocity 
dependences of the parallel and perpendicular components of the diffusion tensor are 

/ 1 v 2 \ ( 3 v 2 \ 

K±(v,r) = (1 - - ^- \Ko(r), K r (v,r) = (1- -—)K (r), (17) 

CSC CSC 

where k is the (isotropic) diffusion coefficient at v — 0. We also justify in the Appendix 
using this velocity dependence as an adequate approximation for the case a = 7/4, for 
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which we find, 

The quantity InAi can be slightly smaller than In A = ln(M cb h/m b h), as discussed in 
Appendix A, depending on the scale over which the diffusion takes place, because scatterings 
with large velocity changes do not contribute to the diffusion rate over a small range of 
velocities. 



3.2. Three Regimes of Capture 

Before identifying the condition for capture from highly eccentric orbits, it will be 
useful to compute the core radius, r c , where the energy-loss time to gravitational radiation 
on a circular orbit is equal to the relaxation time. We define the relaxation time for circular 
orbits as t r>0 = v 2 /[2tr(K,)], where v c = (GM^Jr) 1 / 2 is the circular velocity. The energy 
loss by gravitational radiation for circular orbits is 



d\nE _ 64 G 3 m hh M 2 bh 
dt ~ 5 &r A c 

Equating this to t~Q, we find, 

'GM rhh \ 2/3 1/3 / M chh \ 4/15 



(19) 



- 1 ^7 ( GM <*>A 1/3 ( M cbh \ 

Tc ~ 57 I c 2 ) T ' hh U bh m bh lnAj 



8.6 AU 



M cbh \ 2 / 3 / r bh \ V3 / M cbh 
3 xlO 6 M / V0.7pc/ U.38A" bh m bh lnA/ 1 j 



This radius is extremely small compared to the cluster radius r bh . If black holes were 
captured by diffusing to orbits with a ~ r c , and then losing energy by gravitational waves at 
low eccentricity, the capture rate would therefore also be extremely small. In reality, as we 
shall see now, black holes will be captured from a large range of radii, mostly from highly 
eccentric orbits. We must therefore define a relaxation time scale for these high eccentricity 
orbits. 

To do this, we first calculate the orbit-averaged rate of change of the peribothron 
q = a(l — e) in a highly eccentric orbit. From angular momentum conservation, we have 
v\ = (1 + e)qGM hh /r 2 — > qv 2 sc /r. Hence, 

pd ~W - f d ^r - {^nJ±)°'\, , (2i) 

dt Jo vj 33 \M chh J Vr bh / 
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We then define the eccentric relaxation time i n i(a) = (d(lnq) /dt)' 1 . We now use t r< i to 
demonstrate that there are three regimes of capture: 1) for a < a trans , capture is dominated 
by gravitational radiation, 2) for a tra ns < a < a C rit> it is dominated by direct capture, and 3) 
for a > a CT i t , capture falls off very rapidly and can be ignored. 

A particle orbiting around a Schwarzschild black hole with semi-major axis a much 
greater than the Schwarzschild radius will be directly captured by the black hole if its 
peribothron q, computed by extrapolating the Newtonian orbit, is less than 4 Schwarzschild 
radii (e.g., Misner, Thorne, & Wheeler 1973). When q = 8GM cb h/c 2 , the particle is actually 
brought to 2 Schwarzschild radii by relativistic effects, where the maximum of the effective 
radial potential is located; the particle has then overcome the angular momentum barrier 
and can directly fall into the black hole. Therefore, the phase space density should drop to 
zero below the minimum peribothron 

g mm = ^^0.24AU. (22) 

The black hole can also be captured at larger peribothron if the timescale to change 
the eccentricity by diffusion is longer than the timescale for losing its orbital energy by 
gravitational waves. For l-e<l, the gravitational radiation decay rate is, 

dhiE _ 170 7/2 G 3 m bh M 2 bh 
We determine g m i n , the minimum peribothron that avoids capture, by setting dhxE/dt = 

«-,(«) = °-3H^H l r73Aw^^ J {-) ' (24) 

where the factor 1.73 is the ratio M c bh/(A r bh^bhlnAi) for the values we use in this paper, 
and for InAi = 9.7 (see Appendix A). Hence, the transition from capture by gravitational 
radiation to direct capture occurs at a tran s = 0.35 2 r bh = 17,000 AU. 

Above some critical semimajor axis a cr i t , the diffusion of q over a single period exceeds 
g m i n as given by equation (j22D for direct capture. Hence, during each period, P ~ a 3//2 , 
the black holes are captured with a probability that decreases with semimajor axis as 
(q m i n /a) ~ a -1 , so that the capture rate falls off cx a~ 5//2 . Since the mass within radius 
r increases as r 3 ~ a , captures from orbits with a ^> a cr i t produce a negligible loss of black 
holes. We evaluate a cr i t by setting t r> i = P, and find 

nA , ( N hh yVy M cbh y/ym bh y8/V r bh \^ 
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3.3. Capture Rate Prom The Loss Cylinder 

While k is a function of both v and r, we will solve the diffusion equation at fixed 
r, and temporarily assume that k is independent of v. We will introduce variation in k 
only when we evaluate the loss rate. This is a very good approximation, as we discuss in 
Appendix A. We focus first on the case of a tra ns < a < a cr it, for which capture is typically 
direct as described by equation ([22]) . Since q min is independent of a, the captured orbits at 
fixed r are characterized by a cylinder in velocity space, 



— < = 2^, M < 1, (26) 

^esc * ' ^ ^esc 

where v r and vi are the radial and perpendicular components of the velocity. Note that 
this structure in velocity space is definitely a "loss cylinder" and not a "loss cone" as it is 
often described in the literature. Making use of yet another good approximation described 
in Appendix A, we solve the steady-state diffusion equation fll4l), at fixed v r : 



/(v±;u r ) = fo(v r ) 



l _Hvi/{vi c -vD) 



(27) 



ln(g min /r) 

Hence, the capture rate per unit volume, dC/dV, is given by 

— = - j • dA = 8tt / dv r — , 28 

dv Jcyi-bnd J o m(r/g min ) 

where dA is the area element on the boundary of the capture cylinder ("cyl-bnd"). We 
evaluate this using equations ([15|) and (|17D, 

dC_ 35 (3/4)! In Ax N^ h {Gm hh ) 2 fry 2 _ 3 
dV ~ 16^72 (Vi)! ln(r/g min ) (2GM cbh r bh ) 3 / 2 I W ^ 

= n(r) , 2Q . 

vL ln(r/g min ) ' 

where n{r) = 5N hh / (167irl h ) (r/r bh ) -7 / 4 is the number density of black holes at radius r, 
and 6k,o/v 2 sc is of order the relaxation time. Integrating over volume, the total capture rate 
within some maximum radius r max (determined below) is 

ln(r max /g min ) hh \M chh J r bh P bh 

where k = (2tt)- 1 / 2 (35/8)(3/4)!/(1/4)! ~ 1.76, P bh is the orbital period at a = r bh , and 
where we have made the evaluation treating ln(r/g min ) — > ln(r max /g min ) and lnA x as 
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4/5 

constants. Notice that the capture rate increases as N h ^ as black holes are added to the 
cluster, because r bh oc iV^ 5 , and Pbh oc N®[ 5 . 

As described already by Frank & Rees (1976), the result we have reached in equations 



(f29|) and ( p0|) shows that the capture rate of the black holes is essentially proportional to 
the relaxation time, and the fact that q m \ n / r < 1 increases the time required for the black 
hole to find a capture orbit only logarithmically. The reason is that, over a relaxation 
time, an orbiting black hole will totally change its eccentricity not only due to some large 
deflection in a close encounter, but also due to many small deflections that will change 
the eccentricity by very small amounts, allowing the black hole to effectively sweep over 
all possible eccentricities and find the very narrow range of eccentricity where it can be 
captured. However, this is no longer true for a > a cr i t : when q is brought below q m \ n by the 
random deflections, the black hole will most likely miss being captured unless it happens to 
be at peribothron. 

The maximum radius r max of integration of the capture rate in equation (]30f) must 
therefore be of order a cr i t . For a highly eccentric orbit, the time-averaged radius is 
(r) = (3/2)a. We therefore adopt r max = (3/2)a cri t- Since capture is dominated by black 
holes near r max , we evaluate ln(r/q , min ) there and find, 

ln(r max /g min ) = In ^^T* ^ 13.2. (31) 
16GM cbh 

At the same time, the term InAi can be approximated as (see Appendix A) 

InAi ~lnA--ln ~ 9.7 (32) 

4 8GM cbh 1 ; 



Hence, the ratio of logarithms in equation (p0|) is In Ax/ ln(r max /g min ) ~ 0.73. We are finally 
able to evaluate the capture rate explicitly, 

Since this timescale is much longer than a Hubble time, most of the 24,000 black holes that 
have entered the cluster are still in it and have not been captured. Hence, the actual radius 
of the black hole cluster is close to our initial estimate given by equation (|l"3|). 



Note that we have everywhere used the direct capture formula (22) to calculate q 



rather than the gravitational radiation formula (]24|) , which applies at r < a trans ~ 17, 000 AU. 
Recall, however, that g m i n only enters into the logarithm term. The capture rate from the 
inner part of the cluster is therefore higher than we have assumed, but not dramatically. 
Notice that since q min is a function of a in equation ([24]), the loss structure in velocity space 
is not a cylinder as in the case of direct capture [see eq. Rather, this structure is 

fatter near v r ~ and narrower near v r ~ ±v esc . 
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4. Observable Consequences 

The large number of black holes that move to the center by dynamical friction have 
several observable effects which we now discuss. 

The most important effect is that the stars that formerly resided in the region that 
is now occupied by the cluster of black holes are ejected into orbits at larger radius. 
Therefore, the old stellar population of mass m should have a very large core radius, given 
by r cs = r bh (m bh /m) 2 / 7 ~ 1 — 2pc. This core radius (~ 40") is much larger than the 
value expected from stellar collisions alone, which would only produce a core at the radius 
~ 0.03 pc where the orbital velocities are comparable to the escape velocities from the 
stellar surfaces. The most straightforward test of our model is therefore to measure the 
distribution of low-mass stars and find out if this very large core indeed exists. 

In fact, the bright (K < 15) stars in the inner galaxy exhibit a power-law profile 
(a ~ 7/4) all the way to r ~ 0.1 pc where a core sets in (Genzel et al. 2000; Schmitt 1995). 
However, this is not necessarily inconsistent with a core in the old population: most of the 
observed stars could well be young. Indeed, Genzel et al. (2000) have found that the stars 
they observed very near to the Galactic center are young supergiants, which must have 
formed in a recent (~ lOMyr) burst of star formation. If such bursts of star formation have 
occurred repeatedly, then we would expect that the inner region r < r b h ~ 0.7 pc should 
contain the remains of the bursts from the last t r (r b h) ~ 2 Gyr, but not earlier. Blum, 
Sellgren, & DePoy (1996a) found that the luminosity function in the inner 2' (5 pc) extends 
to several magnitudes brighter than the one in Baade's Window and inferred that young 
stars must be present in significant numbers. Blum, Sellgren, & DePoy (1996b) obtained 
spectra of 19 of these brighter stars. Of those they could age-date, most were younger than 
200 Myr. 

The intermediate-age (~ 2 Gyr) and older stars ejected from the central cluster could 
be expected to be found up to a few r from Sgr A* on "Oort Cloud" like orbits. That 
is, a star from this population would keep getting jolted by the black holes to a more 
and more eccentric orbit until diffusive encounters near apocenter drove it into an orbit 
with a pericenter just beyond r b h. Narayanan, Gould, & Depoy (1996) list 16 presumably 
intermediate-age stars (K < 5) at projected radii of 8 to 20 pc. If these come primarily 
from the central pc, they should be preferentially on radial orbits. Of course, there has 
been recent star formation outside the central pc as well. For example, Blum et al. (1996b) 
estimate the ages IRS 24 and IRS 23 at ~ 100 Myr and ~ 200 Myr respectively, and these 
both lie at a projected distance of 1.7 pc. Hence, not all the Narayanan et al. (1996) stars 
are necessarily ejected. 
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Another interesting consequence of the ejection of low-mass stars from the center is 
that the capture of ordinary stars by supermassive black holes could be much rarer than 
commonly believed. These captures should lead to bright optical flares (e.g., Rees 1988; 
Ulmer 1999), which could be found in supernova searches. A cluster of black holes should 
have formed around all the supermassive black holes in galactic centers by the dynamical 
friction process described here, and these clusters would reduce by a large factor the rate at 
which stars can come close enough to the black hole to be tidally disrupted. Of course, since 
there is a central density cusp containing at least young stars around Sgr A*, there will be 
some tidal captures, but the absence of old stars should reduce the number of expected 
flares. In addition, for galaxies such as ellipticals that are poor in neutral gas, one would 
not expect continuous star formation near the central black holes. Consequently, after the 
black hole cluster had ejected all the old stars, no young stars would form to replace them. 
For these galaxies, flares from stellar captures could be extremely rare. 

We mention also microlensing of a background bulge star as another potentially 
observable effect, although requiring a large improvement in sensitivity and resolution 
of infrared imaging over our current capabilities. If a bulge star at a distance 2kpc 
behind the center could be observed, the angular Einstein radius of Sgr A* would 
be b = 0."8, corresponding to a linear size of te = 0.03 pc. The two images of the 
star could therefore be comfortably resolved. The star would typically take several 
hundred years to complete the "microlensing event" by Sgr A*. These two images would 
then have a microlensing optical depth to be lensed by one of the cluster black holes, 
t ~ (N hh mbh/M B H)(rE/rbh) 5 ^ 4 A ~ 10~ 3 A, where A is the magnification of the image. 
Imaging down to K — 21, one should on average find about one background star at 2kpc 
within an Einstein radius of Sgr A*, and about 100 similar stars near the core radius of 
~ 1 pc (for which the Einstein radius of Sgr A* would be only ~ 0."02). Of course, the 
identification of the two images of a star lensed by Sgr A* would be of enormous interest 
by itself; despite the large number of orbiting black holes in the cluster, the expected rate 
of the "planet-like" events is still low, and several lensed stars would need to be identified 
to have a good chance of detecting the short events over a period of ~ 10 years. 

Work by AG was supported in part by grant AST 97-27520 from the NSF. We are 
happy to acknowledge discussions with Martin J. Rees. 
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A. Diffusion Tensor in a Kepler Potential 

Here, we calculate the diffusion tensor k for the general case of a p oc r~ a density profile 
in a Kepler potential, assuming that the velocity distribution is isotropic. Recall that the 
phase space distribution is given by equations (^) and (|16|). We begin our evaluation at 
v = where the diffusion tensor, Kq, is isotropic. A single encounter with another black 
hole at speed u changes the velocity by Av = 2Gm\^j (bu). Hence the time averaged growth 
of (v 2 ) is 

^-/w*/a(*^)W, 

2(3 - a)a\ N hh (Gm hh ) 2 v 2 csc ( r \V 2 ~ a , A fAn 
~ (a - 1/2)1(1/2)! (2GM cbh r bh ) 3 / 2 WJ 1 ' 1 1 

where Ai is the ratio of maximum to minimum impact parameters, which we evaluate 
below. 

If the diffusion equation is written (for isotropic Kq) as KoV 2 f = df/dt, then 

^ = 6 -V < (A2) 



which combined with equation (|A1|) gives k,q. 

Next, we evalute k±(v) and K r (v) for the special case of a = 3/2 for which the 
phase-space density (eqs. |T3[ and is independent of both radius and velocity and 



consequently the problem is tractable analytically. Up to an overall normalization 7, the 
transverse velocity diffusion is given by 

d(vV) 7 /•"esc r^maxK) /l + cos 2 6<\ (w r + w) 2 rftan 2 6 1 

~ " dtt r / n — i ; — i Tj— , (A3) 



c?t 2 J-v csc J8=o V 2 / |u r + 1>| sec 6* 

where cos 9 ma , x (u r ) = \u r + v \/{v 2 SC + f 2 + 2u r f ) 1//2 . For the total velocity diffusion d{v±)/dt, 
we find a similar expression, but with (1 + cos 2 #)/2 — > 1. We evaluate these expressions 
and find 

d(v 2 ) _ 2 ^i) = 2 a / 

We therefore conclude that k±_{v)/kq = 1 — 0.2(v/v esc ) 2 and K r (v)/Ko = 1 — 0.6(v /v esc ) 2 . 
These results apply to a = 3/2, but we adopt them for a = 7/4 as well. If the true 
coefficient for a = 7/4 is 0.2(1 + e) rather than 0.2, then this introduces an error into the 
capture formula ([30]) of only (l/6)e. Since e is likely to be small, this correction is negligible. 

We now justify two other approximations made in the capture calculation in § 3. First, 
in equation (|27|), we effectively considered the diffusion coefficient as fixed on slices of the 
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velocity sphere perpendicular to the cylinder and passing through v r . In fact, it is fixed on 
spherical shells of constant speed. This "plane parallel" approximation is justified by two 
related considerations. First, k varies very slowing, only by a factor 1.25 over the entire 
velocity sphere. Second, most of the depleted region of velocity space is relatively near the 
cylinder, so that the error in k made by the plane-parallel approximation is extremely small. 

Next, the boundary condition for the solution to the diffusion equation ( p7|) sets the 
phase-space density / at the edge of the velocity sphere (u = v esc ) equal to the unperturbed 
density fo at u = v r . Strictly speaking, if the boundary condition is set at {u = v esc ), then 
one should use fo(v esc ). However, for the a = 7/4 profile, this is zero. (Our approximation 
would be exact for a = 3/2). The basis for our approximation is that the density returns 
essentially to f long before the edge of the velocity sphere, at which point f is not much 
different from fo(v r ). 

Finally, we evaluate Ai = b ma , x /b m i n , the ratio of the maximum to minimum impact 
parameters. This can also be written Ai = e max /e min , where (e min t> esc , e max t> csc ) is the 
range of impulses relevant to the problem at hand. For general relaxation e max = 1 and 
e m in ~ 2Gmbh/(^'Wg SC ) = mbh/M c bh- However, while the harder scatters all contribute 
to the overall relaxation, they do not contribute to diffusion into the loss cone because 
the black hole simply "jumps over" the capture cylinder. Unfortunately, it is not trivial 
to identify exactly what the largest allowed jumps should be: while jumps larger than 
(<?mm/ r ) 1 ^ 2 'Wcsc do not directly lead to capture, they do help maintain the overall velocity 
profile given by equation (|27|) . Note that the profile differs significantly from the background 
for many e-foldings. If the jumps larger than (gmin/^) 1 ^ 2 did not contribute at all, then 
Ai = (<?mm/ r ) 1 ^ 2 A. In view of the distribution's slow approach to the background level near 
the capture cylinder, we adopt A x = A(g min /r) 1//4 = 9.7. 
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